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Abstract 

We extend the maximal unitarity method at two loops to double-box basis integrals with up to 
three external massive legs. We use consistency equations based on the requirement that integrals 
of total derivatives vanish. We obtain unique formulas for the coefficients of the master double-box 
integrals. These formulas can be used either analytically or numerically. 
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I. INTRODUCTION 



A quantitative understanding of backgrounds at the LHC to new-physics signals is im- 
portant to direct searches for physics beyond the Standard Model. A similar understanding 
of signals of heavy Standard- Model particles is important to detailed studies of their proper- 
ties as another window into new physics. Such quantitative understanding requires at least 
next-to-leading order calculations in quantum chromodynamics (QCD). These calculations 
make use of a number of ingredients beyond the tree-level amplitudes required for a leading- 
order (LO) calculation: real-emission corrections, with an additional emitted gluon, or a 
gluon splitting into a quark-ant iquark pair; and virtual one-loop corrections, with a virtual 
gluon or virtual quark in a closed loop. The required one-loop corrections are challenging 
with traditional Feynman-diagram methods, and become considerably more difficult as the 
number of final-state partons (gluons or quarks) grows. 

Some subprocesses are absent at tree level, such as the gluon fusion to diphoton subpro- 
cess, gg — > 77, which is an important background to measurements of the recently-discovered 
new heavy boson (NHB) at the LHC . While these subprocesses are nominally of higher 
order in the strong coupling a s , the large gluon parton density at smaller x can compensate 
for this additional power, giving rise to contributions to cross sections which are comparable 

nn 

to those from tree-level quark- initiated subprocesses pH5(. The production of electroweak 
boson pairs, gg — > Z7, ZZ, W + W~ also falls into this class. Because the basic processes only 
arise from one-loop amplitudes, computation of the NLO corrections requires two-loop am- 
plitudes 6] (as well as one-loop amplitudes with an additional parton, and singular factors 
derived from one-loop amplitudes). 

Two- loop amplitudes are also required for any studies beyond NLO. Next-to- next-leading 
order (NNLO) fixed-order calculations stand at the next frontier of precision QCD calcu- 
lations. The only existing fully-exclusive NNLO jet calculations to date are for three-jet 
production in electron-positron annihilation Q. These are necessary to determine a s to 
1% accuracy from jet data at LEP Q], competitively with other determinations. At the 
LHC, NNLO calculations will be useful for determining an honest theoretical uncertainty 
estimate on NLO calculations, for assessing scale stability in multi-scale processes such as 
W+ multi-jet production, and will also be required for precision measurements of the NHB 
or of new physics once it is discovered. 
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The unitarity method |9H25| . has made many previously-inaccessible one-loop calculations 
feasible. Of particular note are processes with many partons in the final state. The most 
recent development, applying generalized unitarity, allows the method to be applied either 



analytically or purely numerically 



26|-|36| . The numerical formalisms underly recent software 



libraries and programs that have been applied to LHC phenomenology. In this approach, the 
one-loop amplitude in QCD is written as a sum over a set of basis integrals, with coefficients 
that are rational in external spinors, 



Amplitude = coefficient j x Integral^ + Rational . 

jgBasis 



The integral basis for amplitudes with massless internal lines contains box, triangle, and 
bubble integrals in addition to purely rational terms (dropping all terms of 0(e) in the 
dimensional regulator). The coefficients are calculated from products of tree amplitudes, 
typically by performing contour integrals via discrete Fourier projection. 

The unitarity method has also been applied to higher-loop amplitudes. Most prior ap- 
plications made use of 'minimal' generalized unitarity, cutting just enough propagators to 
break apart a higher-loop amplitude into a product of tree amplitudes. Each cut is again a 
product of tree amplitudes, but because not all possible propagators are cut, each general- 
ized cut will correspond to several integrals and their coefficients, and algebra will again be 
required to isolate specific integrals and their coefficients. This approach does not require 
a predetermined general basis of integrals. A number of calculations have been done thi s 
way, primarily in the M = 4 supersymmetric gauge theory and M = 8 supergravity 
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531 ] . but including several four-point ca 



less-than- maximal supersymmetry 54J-I60] 



culations in QCD and supersymmetric theories with 



In maximal generalized unitarity, one cuts as many propagators as possible. In princple, 



this allows one to isolate individua 
of eq. (jl.ip . In a previous paper 



integrals on the right-hand side of the higher-loop analog 



61] , two of the authors showed how to extract the coef- 



ficients of the two master double boxes using a higher-dimensional contour around global 
poles. Each global pole corresponds to cutting all propagators, and addition seeking the 
poles of remaining degrees of freedom (typically generated by the Jacobian from cutting the 
propagators). This approach may be viewed as a generalization to two loops of the work of 
Britto, Cachzo, and Feng [3], and of Forde 20]. In parallel, Mastrolia and Ossola 62] and 



Badger, Frellesvig, and Zhang [63] have showed how to generalize the OPP 18[ approach 



to higher loops. A recent paper by Zhang 64| adds tools required in such an approach for 
reducing integrands to a basis of monomials. 

In this paper, we continue the maximal-unitarity approach of ref. jfllj], taking advantage 
of recent work 65] by Caron-Huot and one of the authors showing that the set of unique 
global poles is smaller than previously assumed. Higher-loop amplitudes can be written 
in a similar form to those at one loop (11. ip . as a sum over an integral basis 66J], along 
with possible rational terms. (At higher loops, the coefficients of the basis integrals are 
no longer rational functions of the external spinors alone, but will depend explicitly on e.) 
Furthermore, at two and higher loops, the number of master integrals for a given topology 
may depend on the number and arrangement of external masses 66]. 

As in ref. 6lj], we use the equations relating generic tensor integrals to basis or master 
integrals in order to ensure the consistency and completeness of the choice of contours. We 
obtain a set of unique projectors or master contours to compute the coefficients of double 
boxes with one or two external masses, one projector for each different master integral. We 
again work only to leading order in e in the coefficients. 

The extraction of the double-box coefficient bears a superficial similarity to the proce- 



dure that would be followed in the leading-singularity approach [45J, |67|, but unlike the 
latter, manifestly ensures the consistency of the extraction with respect to terms that van- 
ish after i nteg ration. Such terms inevitably arise when using the integration-by-parts (IBP) 



approach 68M75I] in relating formally-irreducible tensor integrals to basis integrals. The ex- 
traction of higher-order terms in e or the coefficients of integrals with fewer propagators, 
both of which we leave to future work, would also be different. 

In section [Til we present our parametrization of the loop momenta. In section IIII[ we 
discuss the solutions to the heptacut equations, and in section [IV] we list the global poles 
of the double boxes. We discuss the constraint equations in section [V] and obtain our main 
results for the projectors extracting integral coefficients in section I VI I Section IVHI gives 
examples of applying the technique described here, and provides cross-checks with known 
results. We give our conclusions in section IVIII1 
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II. PARAMETRIZING THE INTEGRAND 



It will be convenient for our pur pos es to adopt a different parametrization for the loop 



momenta from the one used in ref. 61J. The parametrization has the virtue of keeping all 
expressions herein manifestly rational throughout the derivation. In this parametrization, 
we will need spinors corresponding to massive legs. To define appropriate momenta, we 
make use of 'mutually projected' momenta as in the work of Ossola, Papadopoulos and 



Pittau 



1* 



and Forde 



201 ] on the extraction of triangle and bubble coefficients at one loop. 



We use (ki, k 2 ) and (/c 3 , fc 4 ) as mutually-projecting pairs, defining within each pair, 



2/c^i • kj 2 



'i,2 ' 



k 



2 



%2 " %2 2k . 2 . k ^,l- 



By construction, k^ x and M 2 are massless momenta. 
Define 

where i denotes one momentum in the pair, and i the other. 
We note that 

%,1 ' ^j',2 = ' ^i,2 = ' ^j',2 ■ (2-3) 

Define 

7il i 2 = 2*},! • 2 ; (2.4) 

we can solve for it to obtain, 

7,ti2 = %i " %2 ± [(%i • % 2 ) 2 - klklf 2 . (2.5) 

If either momentum kjj is massless, only one solution survives^ ^ j \2 — 1 " 2 • 

With two mutually-projecting pairs, we may have two independent 7s, 712 for the (fci, A; 2 ) 

pair, and 734 for the (fc 3 , fc 4 ) pair. In general, we expect to express our final results in terms 

of the independent invariants Si 2 and S14, the nonzero masses amongst {mi, m 2 , m.3, m.4}, 

along with 7l2 (if mi 7^ 7^ m 2 ), and 734 (if m 3 7^ ^ m 4 ). 

We adopt the following parametrization for the heptacut double-box loop momenta, 

= ^<A 1 -|M|A , r> = i<Ai|/*|A / 1 ] , 

^ = ^(A2|/i|A' 2 ) = ^(A 2 |/i|A' 2 ] . (2.6) 



This form ensures that i\ = = i\. We write the various spinors in terms of the spinors 
corresponding to k\ 2 for £i, and the spinors corresponding to k\^ for £ 2 : 

|Ai + ) = ^' + ) + ^ { ^\2^) , 

[4> 1 b l 

IV) = 6-^|||3 b ' + )+e 4 |4 b ' + ), (2.7) 
[1 b 4>1 

|A' 2 -> = ^| P 3fj|3 b '->+^|4 b '->. 

Without loss of generality, we can set £i = 1 = £4. We will use the simplified notation, 
\j] = and |j) = in later sections of the paper. 

III. HEPTACUT EQUATIONS 




k\ hi 



FIG. 1: The double box integral P 2 * 2 . 

We focus on the planar double-box integral, shown in fig. [TJ which is given by, 

p „_ fd^d^ 1 

2 < 2 J (2n) D (2n) D - k x )\i x - K X2 )\t x + £ 2 ) 2 ^ _ fc 4 )2( 4 _ ^ 34 ) 2 ' 1 ' ] 



where Ki...j = hi + ■ ■ - + kj, and the notation follows ref. 66] . We allow any two of the 
first three external legs to be massive, leaving the three- and four-mass cases for future 
presentation. We leave the Feynman is prescription implicit. 

The heptacut equations correspond to cutting all seven propagators, that is setting all 
denominator factors to zero simultaneously. Imposing the heptacut can be implemented by 



replacing the original contour of integration for the loop momenta, running along the real 
axis for each component of £^ 2 , by a seven-torus encircling the solution surface(s) for the 
vanishing of these denominator factors. Two of these equations are already solved by the 
parametrization in eq. (I2.6p . Using that parametrization, the remaining heptacut equations 
involving only £\ take the form, 

(d - k{f = =* 

(h - K a f = =* (3.2) 

= (fi+ftjj&fi) 712 ' 
In the general case {k\ ^ 7^ these equations have the solution, 

pl 712^12 - (712 + kpk\ = -, 

^1 ^2 _ pL2 — SI 7 

/12 ^1^2 

fc 2 ( gl2 - 7l2 -A; 2 )4-^ = g 

^ 2 " 6(7i 2 2 - • *i ~ 6 ' 1 ' ] 

There are similar solutions for the equations involving only £ 2 . 

The first of these solutions has the correct limits when either k\ — > or k\ — > (or both); 
the second solution, in contrast, is replaced by the pair of solutions, 

6 = 0, £ 2 free; G r £ 2 = , 6 free . (3.4) 

More explicitly, when k\ = but k\ ^ 0, 

^ = — = l + p 2 ; (3.5) 

712 

when fc| = (whether k\ vanishes or not), 

£ = !■ (3-6) 



The last heptacut equation is £\ ■ £2 = 0, which takes the form, 

(l b 4 b ) (2 b 3 b ) 
(2 b 4 b )(l b 3 b > 



= (l b 4 b )[l b 4 b ] ^1 + 6 + 6 + 66 



V^ 4+ 6 £3 6e 3 [2 b 4 b ][l b 3 b lJ- (J - 7j 
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This equation has two solutions corresponding to the two factors, 

1 + 6 



(3.8) 



where 



r 



In fact, 



r = r 



(l b 4 b ) (2 b 3 b ) 
(2 b 4 b )(l b 3 b > 
[l b 4 b ] [2 b 3 b ] 
[2 b 4 b ] [l b 3 b ] ' 



734(712 + ml) 



(3.9) 



(734 + ml) [(712734 - mfm|) (712734 - m\m\) + 712734S12S14] 
x [(712 + ml) (734 + m^wij +mg- s 14 ) + 2777,2/77,3812 

+ (712 + ml)ml(m\ -m 2 2 - s 12 ) + (734 + ml)m 2 2 (m\ -m\ - s 12 )} 



(3.10) 



As in ref. 65], we must consider three different classes of solutions, corresponding to 
different configurations of the external masses. 

2 3 



FIG. 2: The short-side two-mass double-box. Single lines indicate massless legs, and doubled lines 
massive legs (here the two legs on the left-hand side). 

In considering massless limits, there are three classes of cases to consider: (a) when neither 
pair (1,2) or (3,4) contains a massless momentum; (b) when exactly one pair contains 
a massless momentum; (c) when both pairs contain massless momenta. The first class 
corresponds to the four-mass double box; the second class, to the three-mass double box and 
the short-side two-mass one, shown in fig. EJ the last class, to the remaining configurations 
- diagonal and long-side two-mass double boxes, shown in fig. 131 along with the one-mass 
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^ 2 









^ 2 







(a) 



(b) 



FIG. 3: Other two-mass double-boxes: (a) the diagonal one (b) the long-side one. Single lines 
indicate massless legs, and doubled lines massive legs (in (a), the lower-left and upper- right legs, 
and in (b), the two lower legs). 

and massless ones. Case (a) is given by the solutions above, but we will not consider it any 
further in the present paper. 

In class (b), it suffices to take the massless momentum to be either k 3 or k 4 (or both). 
There are two solutions to the £ 2 equations to consider, £3 = and £3 = 0. For the first, the 
last heptacut equation becomes, 



= (l b 4 b > [l b 4 b ] (1 + 6) U'iC + ^ + m + r-fl) 

S"2 C9 / 



-£2^3 

6 



(3.11) 



This equation again has two solutions, £ 2 = _ 1 (with £ 3 being identified as the remaining 
variable z), and the other having, 



£3 = 



(3.12) 



£6 + fi' 2 

with £ 2 = z. The second solution to the £ 2 heptacut equations has £3 = 0, whereupon the 
last heptacut equation becomes, 



= <l b 4 b > [l b 4 b ] (1 + ^ + & + rfc&) + 
which also has two solutions, 

6 = - >, 
si 

1 + 6 

l + r6 



6 



3 



(with £3 — 2)) and 



(3.13) 



(3.14) 



(3.15) 
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(with £ 2 = z). Overall in class (b), we have four solutions to the heptacut equations, in 



agreement with the discussion in ref. 



65|, with = £[, £ 4 = £ 4 in all, and the four variables 



(£2, £2; £3; £3) taking on the following values, 



S 1 : 


' £2 


S 2 : 


(4 


S3: (-1,- 


<S 4 : 


(4 



1 + ^,o 



1 + TZ' 



(3.16) 



£> + t£ 2 

In class (c), we have four solutions to the pure-£i j2 equations to consider, 



6 = and £3 = 

6 = and £ 3 = 0. 

£ 2 = and 6 = 0, (3.17) 

£ = and £ 3 = 0. 

In the first of these solutions, the last heptacut equation simplifies to, 

= <l b 4 b > [l b 4 b ] + £ 2 £ 4 + &s + r£ 2 £ 3 ) , (3-18) 

which has the solution, 

t> _ £U£l + £2) /o 1QN 

(with £ 2 = 2). In the second solution in eq. f)3.17p . the last heptacut equation simplifies to, 

0=(l b 4 b )[l b 4 b ] (1 + 6)^4 + ^4) (3.20) 

which has two solutions, 6 = — 1 (with £ 2 = z), and £ 2 = — £j (with 6 = 2). 
In the third solution in eq. (I3.17p . the last heptacut equation simplifies to, 

= (l b 4 b >[l b 4 b ](l + 6) 

x(&4 + &'s), (3.21) 

which also has two solutions, 6 = — 1 (with £3 = 2), and £3 = — £ 4 (with £ 2 = 2). 
In the last solution to eq. (I3.17p . the last heptacut equation simplifies to, 

= (l b 4 b ) [l b 4 b ] ££(1 + 6 + £3 + r£ 2 £ 3 ) , (3.22) 
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which has a single solution, 
(with ^2 = z). 



Ovcra 
refs. 61 



1 in class (c), we have six solutions to the heptacut equations, in agreement with 



65[, with £i = £[, £ 4 = £ 4 in all, and the four variables (£2, £2; 63> £3) taking on the 



following values, 



5 i: (0,-^,2,0), 

S 2 : (0,2,-1,0), 

S 3 : (-1,0,0,^), 

S 4 : (2,0,0,-^4), (3.24) 

5fi : (<>,*, o, 



IV. GLOBAL POLES 

As we take the cuts strictly in four dimensions, the two loop momenta have a total of eight 
degrees of freedom. After changing variables to the various £j and replacing the original 
contours of integration with tori encircling the locations of the zeros of the denominator of 
the integrand localizes it on the solutions of heptacut equations. We found the solutions to 
these equations in the previous section. In the all-massless, one-mass, diagonal two-mass, 
and long-side two-mass double boxes, corresponding to class (c), there are six independent 
solutions. In the short-side two- mass double box, corresponding to class (b), there are four 
independent solutions. 

In either case, the heptacut equations fix only seven of the eight degrees of freedom. 
Performing the corresponding contour integrals gives rise to an inverse Jacobian from di- 
agonalizing and linearizing the denominator in the neighborhood of the simultaneous zeros 
in each of the denominator factors. In all cases, this inverse Jacobian also has poles in the 
remaining degree of freedom z. In certain solutions, powers of one or the other of the loop 
momenta will also have poles in z. These poles make it possible to choose contours for the 
z integration which further localizes the integral. As we shall see, this in turn makes it 
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possible to distinguish between the various master integrals. 

There are two distinct Jacobians that arise. One emerges from the change of variables 
ii2 t° the £ parameters. This Jacobian appears in the numerator. The other is the inverse 
Jacobian from performing the heptacut contour integrals. Because the parametrization in 
the previous section automatically ensures that £f = = £ 2 , it is not suitable for computing 
these Jacobians. As a remedy, add an additional vector to the original parametrization (12.61) . 

t x = (ArHA'r) + a?7i, 

% = (A 2 -|/i|A' 2 -) + C 2 %, 

where Q are complex numbers, and the r]i are null vectors satisfying ^ 1 |A 1 + ) / 7^ ^JA'i - ) 
and ^|A 2 +)^0^ 2 |A' 2 -). 

In both class (b) and (c), we may choose 771 = k 2 and 772 = k\. This is a suitable choice 
so long as 7^ 0, so that the non-collinearity condition above is satisfied. For generic 
momenta, this holds. With the additional term, the equations l\ = = £ 2 become, 

= 7i2^Ci , 

= 734^2- (4.2) 

The requirement that ^ ^4 leaves only the obvious solutions, Ci = and (2 = 0, 
respectively, to these equations. 

The Jacobian for the change of variables is, 

dt 

Jj = det — ^- , (4.3) 

where v jA = Q, v h2 = £[, v %2 = v i,3 = £2, ^2,3 = £3, v 1A = £ 2 , and v 2A = £ 3 . 
In all solutions in both classes (b) and (c), these Jacobians have the values, 

7l2^1 ' ^4^1 



Jl 



The Jacobian from performing the heptacut contour integrals (which will appear in the 
denominator) is, 

BP 

J h = det |-L , (4.5) 
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where Pj is the j-th propagator denominator, and the set {wj} are the seven variables (out 
of the eight £[, £2, ^ £3, £3; £i, Cii C2) frozen by the contour integrals. 

This Jacobian is not the same for all solutions. As discussed earlier, it will appear in the 
denominator; the overall Jacobian appearing after performing the heptacut contour integrals 
is, 

J? = 4^ . (4.6) 



For our purposes, it will be convenient to factor out an overall prefactor which would oth- 
erwise cancel in final formulae for integral coefficients, defining the reduced Jacobian for 
class (b) as, 



and for class (c) as, 



where 



Js, = J\ 



(m £2)^4 



Jf, = J 



7* 



£1^4 



7* 



J 7* 

712734 



34 



(4.7) 



(4.8) 



(4.9) 



OiOOO 




.Sq 000 



omo 



FIG. 4: A representation of the solution space for the class (b) heptacut equations, showing the 
four independent solutions Si, and the locations of the eight global poles Qi. 
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Within class (b), we then find, 



J a, 



(6,1) 



(£l - 



(6,2) 



2 



J 



(6,3) 



+ 1)£ + £ 2 ) ' 



The first of these Jacobians has poles at, 



the second, at 



and the third, at 



z 







z 
z 



•i, 

a. 



solution «Si , 
solutions <S 2 ,4 , 
solution S3 . 



(4.10) 



(4.11) 



(4.12) 







£1 £2 



(4.13) 



£i - r£ 

Integrals with powers of the loop momentum in the numerator will have powers of the 
various coefficients £2, £25 £3, £3 once re-expressed in terms of our parametrization of loop 
momenta ( 12. 6p . Some of these coefficients also have poles in z, which leads to additional 
poles in the integrand beyond those in the Jacobian above. From eq. f!3.16p . we see that the 
additional poles are at, 



0. 



T 



(4.14) 



for solution S 2 , and, 



0. 



£i 



(4.15) 



for solution Sa 
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Poles in z complete the poles in the other seven degrees of freedom to provide a global 
pole. One might expect each pole in z within each solution to yield an distinct global pole; 



but this is not the case. As pointed out [65] by Caron-Huot and one of the authors, global 
poles within one solution may in fact end up being identical to global poles in another, in 
the sense of being located at the same values of After removing such duplicates, one 
finds a total of eight global poles. At all eight, = £[, £ 4 = £ 4 , and 0=0 = (2, while the 
values of (£2, £2' £3' £3) are as follows, 

r> ( £2 s (2^(1 n 

61 ■ ri'^'i^S' 



ft = (-jU.0.0) 



£3: 

G A : (-1,-^,0,0), (4.16) 



5 




Q 7 : f--j-,-"r£ 2 ,oo,0^ 
G 8 : (0,oo,-l,0). 



We obtain the same poles whether we pick 7 12 = 7^ or 7i2> so we can choose either. In 
future numerical applications, it may be more appropriate to average over the two choices. 

In the poles where one parameter or another goes to infinity, we must specify the order 
of limits carefully. 

In global pole Q 3 , we should set, 

£2 = ^ , £3 = , & = - , 6 = s , * = , (4.17) 

z r£ 2 + z£[ 

followed by the limit z — > 0; in global pole Q$, we should set, 

£ = ^, 6 = 0, & = Ji -f[ + P , 6 = z, (4.18) 
followed by the limit z — > — f£ 2 /£[] in global pole Qj, we should set, 

& = ^,& = 0,&= — ^Z7. 6 = z, (4.19) 

Z TZ + 1 
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followed by the limit z — > 0; and in global pole G$, we should set, 



(4.20) 



z'™ ^ rz + l 
followed by the limit z = —1/r. The solutions and global poles are depicted graphically in 

fig.H 



6/ / J\ 


G 3 




uQi . / 
1 1 . / 







oo< 



5s •oo 



FIG. 5: A representation of the solution space for the class (c) heptacut equations, showing the 
six independent solutions Si, and the locations of the eight global poles Q{. 



Within class (c), we find, 



J 



M) 



j 



(c,2) 



J 



(c,3) 



solutions Si 2 , 



solution ^4 



Solutions 5 6 • 



The first of these Jacobians has poles at, 



the second, at 



z 
z 

z 
z 



-3; 



(4.21) 



(4.22) 



(4.23) 
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and the third, at 

z = 0, 

z = -1. (4.24) 
Here, looking at eq. (I3.24p . we see that powers of the loop momentum in the numerator 



will induce poles at, 

1 



T 



(4.25) 



for solution <S 5 , and, 



z = (4.26) 

T 



for solution <S 6 . 

Once again, once the dust settles on comparing ostensibly different solutions, one finds 
eight independent global poles, with 3 = 3, 3 = 3' Ci = = C2, and the following values 
of (6,3,6,3), 

(-1,0,0,0), 
& : (0,-£,-l,0), 
£ 3 : (0,-^,0,0), 

£4: (0,0,0,-3), (4-27) 

£ 5 : (-1,0,0,-3), 

£ 6 : (o,-|,0,oo 

£ 7 : f-~,0, oo,0 
£ 8 : (0,0,-1,0). 

Once again, for the global poles with an infinite coefficient, we must take the limit care- 
fully; for global pole Q 3 , we should set, 

3 = 0,3 = 0,6 = —^7,6 = z, (4.28) 
and then take the limit z — > — 1/r, while for global pole C/s, we should set, 

6 = 0,6 = 0,3 = -1|±^ , 3 = z, (4.29) 

rz + 3 

and then take the limit z —> —£[/f. The solutions and global poles are depicted graphically 
in fig. [5j 
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V. CONSTRAINT EQUATIONS 



As we have seen in the previous section, each double-box integral with a general numerator 
has eight independent global poles. One might therefore imagine obtaining eight independent 
equations for coefficients of integrals by integrating both the left- and right-hand sides of 
the two-loop version of eq. (11. ip on a T 8 contour surrounding each global pole in turn, 

/ ^Amplitude Integrand = V Cj i d % ah ^Integrand, . (5.1) 

where T 8 (Qi) denotes an eight-torus encircling the ith global pole, and where is the 
corresponding Jacobian derived in the previous section. The v a j,, as before, denote the eight 
variables 0,2, 6,3, and £1,2,3,4- 

These equations would not all be valid, however, because the derivation of the reduced 
form on the right-hand side assumes that a number of integrals, with certain kinds of nu- 
merators, vanish identically. We must restrict ourselves to contours, whether encircling one 
global pole or multiple ones, that enforce the vanishing of these integrals. As illustrated 



in detail in a simple one-loop example in ref. [61], non- vanishing numerators which yield 

a vanishing integral over the standard real contour will generically not vanish on a linear 

combination of contours encircling global poles. 

We focus on the double-box topology. Denote the insertion of the polynomial f{£i,£ 2 ) 

into the numerator of the double box by, 

P ** lf(f n]= fd^d^ f(£i,£2) 

2, 2 LM i, 2jj j (2 7r )^( 27r )^^(£ 1 -A: 1 ) 2 (^i-^i2) 2 (4+4) 2 ^(4-A:4) 2 (^-^34) 2 ' 

(5.2) 

For class (c), there are two master integrals, which we take to be 

I t = p**[l] , and I 2 = P;* 2 [U ■ k 4 ] . (5.3) 

For class (b), there are three master integrals, which we take to be, 

h = P% [1] , h = P% [4 • h] , and h = Pll [£2 ■ h] . (5.4) 

Vanishing integrals with non-vanishing numerators can all be thought of as integrals of 
total derivatives. They fall into two distinct classes, parity-odd and parity-even. The former 
consist of five integrals with insertions of Levi-Civita symbols, 

P 2 *; 2 [e(€i, k 2 , h, h)) , P%[e{l 2 , k 2 , k 3 , h)} , P^Wu £ 2 , h, k 2 )) , 
P;; 2 [e(£u£2,k 1 ,k 3 )}, and P%[e{h, £ 2 , k 2 , k 3 )) . (5.5) 
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whose vanishing on our chosen contours imposes one set of constraints. 

The parity-even total derivatives are precisely the integrands which yield the equations 
required to fully reduce all possible integrals arising in a gauge theory to the two or three 
master integrals. 

Imagine computing a two-loop amplitude with four external momenta diagrammatically. 
At the first step, the numerator of the two-loop double box will contain dot products of the 
two loop momenta £\^ 2 with external polarization vectors, external momenta, and external 
spinor strings (if any external legs are fermions). We can express all dot products of loop 
momenta with external vectors in terms of eight dot products: £j ■ ki, £j ■ k 2 , £j ■ k^, and £j ■ v, 
where v M = k%, k 2 , k^). Just as at one loop, odd powers of v will give rise to vanishing 
integrals, as expressed in the Levi-Civita constraints discussed above. Even powers can 
again be re-expressed in terms of the other dot products (up to terms involving the (— 2e)- 
dimensional components of the loop momentum). All integrals can then be rewritten in 
terms of the six dot products of the loop momenta with the external momenta. 

Of these six dot products, three of them — £\ ■ k\, t\ ■ k 2 , ti ■ k± — can be rewritten as 
linear combinations of the propagator denominators and external invariants. One additional 
dot product of £ 2 — say £ 2 -k 2 — can be rewritten in terms of the remaining two (£i ■ k^ and 
£ 2 - ki), propagator denominators, and external invariants. The remaining two dot products 
are called irreducible. At a first stage, then, before using IBP identities, we can reduce an 
arbitrary double-box integral appearing in a gauge-theory amplitude to a linear combination 
of the 22 different integrals that can arise with powers of the two irreducible numerators. 

At first glance, all integrals of the form, 

P** 2 [(£ 2 ■ h) m (h-h) n ], (5.6) 

would appear to be irreducible. In a general gauge theory, all integrals with < m + n < 
6, m, n < 4 can appear, giving a total of 22 integrals. The IBP relations will allow us to 
reduce the integrals to a small set of master integrals. The number of IBP relations depends 
on the configuration of external momenta. For class (b), we have 19 such relations, leaving 
the three masters in eq. ( 15.4[) . For class (c), we have 20 such relations, leaving us with the 
two masters in eq. (15.31) . Each IBP relation gives us a nontrivial numerator which integrates 



to zero; explicit examples may be found in ref. 



fill 



In order to ensure that these equations continue to hold after performing the global 
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contour integrals, we seek a linear combination of T 8 contours surrounding the different 
global poles, 



d 8 v Jf<F , 



(5.7) 



and replace F in turn by each of the parity-odd or -even numerators discussed above. This 
gives us five equations for the parity-odd numerators, and either 19 or 20 for the parity-even 
ones. 

It turns out that one of the five Levi-Civita equations is linearly dependent on the others, 
just as in the purely massless case considered in ref. {fy]]. This leaves us with four independent 
constraints; in class (b), these have the form, 



/ -1 1 -1 o\ 
0-1100 
0-1011 
\-l 1 1 000/ 



and in class (c), 



/ -1 1 oo\ 
1 -1 
-110 
\ 00 -11/ 



a 3 
<%i 
a 5 
a 6 

a 2 

Cl{ 
«5 

a 7 



0; 



0; 



(5.8) 



(5.9) 
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11111-10 



0; 



1 1 1 1 
110 11-1-1 



0. 



(5.10) 



Similarly, only one IBP equation is independent for class (b), 

/ai\ 

a 3 

a .4 

«5 
«6 

and only two for class (c), 

/ ax\ 

a 2 

a i 

a 5 
a 6 
a 7 

The simultaneous solutions to these equations are choices of contours suitable for inte- 
grating on both sides of eq. (11. ip in order to derive useful relations for coefficients. 



(5.11) 



VI. PROJECTORS 



Imposing the constraints obtained in the previous section will leave us with three non- 
trivial linearly independent combinations of contours for external kinematics in class (b), 
and two combinations in class (c). Using these contours on the two- loop version of eq. (II. ip 
restricted to four-point amplitudes will then give us three or two equations, respectively, 
for the coefficients of the three or two different master integrals (15. 4p and (15. 3p . We can 
diagonalize these equations by requiring in turn the vanishing of each master integral along 
the contour. The resulting equations for the coefficients take the form, 

8 . 6 

a = J2 a if s d%v ^ J §d Yl II a p ] M , (6-1) 

j = l JT 8 (Gj) particles p=l 

helicities 
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where the product in the integrand is over the tree amplitudes at all vertices, with their 
dependence coming through the parametrization (12. 6p . The product is then summed over 
all allowed choices of internal particles and helicities. 

This abstract form of a multidimensional contour integration hides a subtlety whose full 
exploration we postpone to a future paper. In multivariate complex contour integration, 
the result is not necessarily independent of the order of integration. In our context, this is 
connected with accidental sharing of singularities between different integrals, a feature not 
present in the one-loop case. Once we specify the ordering, the accidental sharing disappears, 
as does any ambiguity in the result. In the present context, the order in eq. (16. ip is specified 
by performing the z integration last, 

(6.2) 



d v J = <h dz JfjZ 



'T s (gj) Jc s (z°) 
where the integrand is first to be evaluated a heptacut solution S(Qj) containing the global 
pole, and the contour integral in z is then to be performed on a small circle C$ surrounding 
the value of z at the given global pole. (In general, there is more than one heptacut solution 
containing a given global pole; we can pick either of the two solutions, so long as we pick 
the corresponding Jacobian Jfj-) 

Diagonalizing the equations, we obtain projectors or master contours isolating a single 
master integral. Each is specified by a vector giving the values of the a« to be used in the 
integrations above. In class (b), the projectors for the three masters are, 

/.= Pi=; (-1,-1, o,o, 2 i { f - sa) - 1. 2 i ( r " ml) + 1) ■ 

N (b) 

h- P 2 = ^- (1,-1, -1,1, 2, 2, 1,1), (6.3) 
N (b) 

I 3 : P 3 = ^(0, 0,0, 0,0, 0,-1,1) , 



where, 



T {b) 2( 7l 2 2 - m\ml) 



/V 1 

2 — / 2 2\ ' 

7i2(si2Si4 -m{m%) 

iVf = - A ^- m p 2 . (6 .4) 
S12S14 — m^mg 

The formulae in previous sections allow for any pattern of masses within class (b); but the 
projectors here are given for the special case that = 0, ni\ 7^ 7^ 777.2, and hold whether 
m 3 is zero or not. 
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In class (c), they are 

t D _ ^i ( 1 -. i 1 2m 4 ( s i2 - m?) 2m^ (sis - m 2 

Ji : -Ti — — — —1, —1, 1, 1 



4 \ ' S12 (su - m|) + mfm|' s 12 (s i4 - mf) + m\m\' 

2m\ (s 12 - mf) 2m 2 (s 12 - mf ) 



(6.5) 



/ 2 : P 2 = ^- (-1,-1, 1,1, -2, -2, -3, -3) , 
where 

S12S14 

m c} = — - — (1+ 4 1 — — ). 6.6 

As above, the formulae in previous sections allow for any pattern of masses within class (c); 
but the projectors here are given for the special case that m 2 = 0, mi ^ 0, and at most one 
of m 3 and m 4 is non-zero. 



VII. EXAMPLES: HIGGS AMPLITUDES AND FORM FACTORS 

In this section we will work out cut expressions, residues and integral coefficients for the 
one-mass double box integral in a few simple but nontrivial four-point examples. We use 
these examples to illustrate how to apply the method described in previous sections. In 
particular we discuss a single massive scalar amplitude and a form factor. A simple way to 
introduce a massive scalar in a gauge invariant way is to add the following interaction to 
the pure Yang-Mills theory: 

Ant. = -^Tr(F^), (7.1) 

where A denotes the scale at which the physics underlying this higher- dimension operator 
has been integrated out, and if is a massive scalar boson. This is the simplest gluon-scalar 
interaction for a colorless scalar not carrying a conserved charge. It is precisely the form 

n 

of the effective Lagrangian [76j for the Standard Model Higgs-gluon interaction, once the 
top quark has been integrated out. Recently, Gehrmann et al. have computed the two-loop 



H — > ggg amplitude in this effective theory [77J] , and simplified by Duhr [78| 
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The computation of the scalar's three-gluon decay amplitude can also be viewed as the 
computation of the three-gluon form factor for the operator, 

C» g = -lTr(F^). (7.2) 

We can compute such form factors using the same approach as for on-shell amplitudes, 

(0\O(q)\l,2,...,n) = A°(q,l,2,...,n). (7.3) 

In addition to the operator O g , we will also consider a well-studied scalar operator in J\f = 4 
super- Yang-Mills theory (sYM), 

a^Tr(0 12 12 ), (7.4) 

where the subscripts refer to .R-charge indices. These indices will not play a role here; one 
can think of 0i 2 as simply a complex scalar with a conserved U(l) charge: 12 = + = 
012 = 0~ = 4>. Brandhuber, Travaglini and Yang recently computed the three-particle form 



factor of this operator at two loops [79|. We compare our calculation of the double box 
coefficients to their results. 

For our computations, the only tree ingredients we need are the three-point amplitudes 
and two-particle form factors. The nonvanishing form factors, after stripping off color factors, 
are, 

(1+ 2+|CgO) = [1 2] 2 , (I" 2-\O g \0) = (1 2) 2 , (1+ 2+|a|0) = 1 . (7.5) 
(We omit an overall dimensionful factor of 1/A.) 



A. Computing the integral coefficients 

Let us follow the procedure derived in previous sections in several steps. First we form the 
product of amplitudes appearing in eq. (16. ip . evaluated at each of the heptacut solutions 
SiSe, depicted in fig. |6j This is the integrand of our remaining contour integral. This 
corresponds to evaluating the factor I in eq. (16. 2 p at each of the solutions S where it will 
be needed. 

Next, we evaluate the remaining contour integral in z by residues. Here, we should 
distinguish between two types of poles: those arising from the heptacut Jacobian (l4.22fl4T24|) : 
and those arising from powers of the loop momentum in the integrand itself (I4.25fj4~26]) . The 
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I 






FIG. 6: The six solutions to the heptacut kinematics of a one-mass double box. White and black 
blobs indicate chiral and antichiral kinematics at the vertices, respectively. The gray blobs denote 
vertices with no definite chirality. 

integrand can safely be evaluated at the former type, where it is finite. The overall residues 
are easily obtained by multiplying the integrand by the residues of the Jacobian, which only 
has simple poles. 

The poles arising from powers of the loop momentum are at infinite values of one of 
the loop momenta. (The poles are nonetheless at finite values of the remaining integration 
variable z.) The Jacobian factor remains finite at these poles. Here the procedure is more 
involved, because in general the poles may be of higher order. The order of the pole depends 
on the ultraviolet power-counting behavior of the process in the given theory. In M = 4 
sYM, and for the form factors we consider here, the poles will be simple. A general fermion 
and scalar content in the theory would lead to poles of higher order. (From a computational 
point of view, the residues at infinite loop momentum require a more careful treatment 
as they may receive contributions that arise from the interference between singular and 
finite terms in the Laurent expansion of the various loop-momentum factors entering the 
integrand.) 

For each kinematical solution Sj we denote the integrand by Xj, 

6 



particles %=\ 
hclicitics 



(7.6) 



At each global pole, the remaining integrand in z has a pole at a single location z®\ we 
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evaluate the residue at that pole, 

Rg. = ResJ /)i J(e i ), (7.7) 

as described above. 

The results below are valid for any massless gauge theory with nj fermions and n s (com- 
plex) scalars, both in the adjoint representation. The kinematical solutions and poles in 
the one-mass case are illustrated in fig. [6] and fig. El respectively. We can read off the 
amplitudes from the pattern of white and black blobs in the figures. Maximally helicity vi- 
olating (MHV) amplitudes are non- vanishing for chiral kinematics, but vanish for antichiral 
kinematics. The white blobs therefore give rise to MHV amplitudes. Similarly, MHV am- 
plitudes are non-vanishing for antichiral kinematics, but vanish for chiral kinematics. The 
black blobs therefore give rise to MHV amplitudes. More detailed information on notation 
and the location of the poles in the one-mass case is summarized in table HI 



2 3 2 3 2 3 




2 3 2 3 2 3 




FIG. 7: The Jacobian poles are located at the intersection points of the six kinematical solutions. 
They correspond to situations where vertical rungs becomes soft (in the massless case), as indicated 
by the absent lines, or in a more general situation where kinematics becomes collinear, as indicated 
by the dashed blobs. 

1. The (O g , 2~ , 3~ , 4~) Configuration 



It will be helpful to label the internal cut lines as follows, 

4 = 4 - h , 4 = 4 - h - k 2 , 4 = 4 - h , 4 = 4 - h - h , l 7 = -4 - 4 • (7.8) 
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TABLE I: The properties of the poles in the one-mass case: Here we assume that k\ = fcf = k\ 



0. 



The notation inj indicates that the pole is located on the intersection point of the Riemann spheres 
labeled by Si and Sj, and com corresponds to poles in Si at infinite 4 momenta (see [(35]). 



pole 


location 




4 


4 


Gi 


1 n2 


(0,-1,-1,0) 




I[4l711 b )#f 


G2 


3n4 


(-1,0,0,-1) 




l<4[7^]fl 


G 3 


6ni 


(0,-1,0,0) 






Ga 


5n3 


(-1,0,0,0) 






G 5 


2n5 


(0,0,-1,0) 


k\ 




Ge 


4n6 


(0,0,0,-1) 


k\ 




G 7 


oo 6R 


/ n [31 b ][4 2l n \ 

(°> Mil*]* ' 00 ) 




oo[3|7^|4) 


Gs 


oo 5R 


\ (32)(41 b )' U '°°' U ^ 


-|(3|7^|l b ]#| 


oo(3|7^|4] 



For each heptacut solution, compute the cut integrand by multiplying together tree ampli- 
tudes for each vertex, and summing over different possible assignments of internal states. 
For example, for the evaluation on Sq, here only a single configuration of internal helicity 
assignments contributes, 

x 6 = H)M Hi ^ ( iH A > -ti) A v>(2-,ez, -tf M<°>(3- -i-,i+)A(°\±-, -it, it) 

xAV>(et,£f,£t)AW(-ei, -£e, -Q) 

i2 (24) 3 (3 4) 3 [hhf (4 4) 3 [44] 3 



1 13 



y 



(44) (42) (44) (43) [44] [44] (44) (44) [44] [V 



7j 



,4<2| 14 14) (43) 



= < [23] L > (7 - 9) 

where the last line is obtained after using some spinor manipulations. 

The integrand takes the same simple form on <Si and £2, whereas it vanishes identically 
on £5, <S>3 and S4, that is 

4 (2|14|4) (43) 

X 6 =Xi=X 2 = m 1 — , X 5 =X 3 =X 4 = 0. (7.10) 

These results are the same independent of the number of scalars or fermions in the gauge 
theory. In particular this means that, for the (O g , 2", 3~, 4~) configuration, all (massless) 
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gauge theories have the same double-box cuts as the M = 4 theory, explaining the simple 
dependence on loop momenta in eq. (17.101) . 

To obtain the eight residues we should evaluate the integrand in the neighborhood of the 
global poles, obtaining 



^6,1(^3) = ^1,2(^1) =^lS 3 4 



4. (23) (34) 



[2 4] ' 
4 (2 4) (4 3) 

= ^1*34 J23] ' 

X 5 (£ 8 ) = X 2 , 5 (£ 5 ) = Z 5i3 (£ 4 ) = X 3 ,4(£ 2 ) = X 4 , 6 (^ 6 ) = . (7.11) 

where the notation X it j indicates that either solution Si or solution Sj can be used. The 
integrand happens to be finite at all global poles, so only the Jacobian poles can give rise 
to nonvanishing residues. As mentioned, the Jacobian has only simple poles, ~ 1/z, at 
the intersecting points of the kinematical solutions. With our normalization, there is no 
additional dimensionful quantity that appears (it was removed in the transition from to 
J § ). 

Finally, after multiplying by the integrand points Xj by the Jacobian factors there are 
exactly two nonzero residues: 

Rg 3 = -Rg, = . (7.12) 

Note that Rg 7 = even though the integrand is nonvanishing at this point. In vector 
notation the residues are 

R = m i g 34 ^] (-1, 0, 1, 0, 0, 0, 0, 0) . (7.13) 

The integral coefficients are then given by a sum over appropriately weighted residues. 
For the double box with unit numerator we have 

ci = Pi • R = -^{-Rg 1 - Rg 2 + Rg 3 + Rg 4 ~ Rg? ~ Rq s ) = -^.§34 ^ ^ - , (7.14) 
and for the double box with loop-momentum numerator 

c 2 = P 2 ■ R = = - l , (-Rg 1 - Rg 2 + Rg 3 + Rg± ~ 2Rg 5 ~ 2Rg a - 3Rg T - 3i?g 8 ) 

- l m *a (23) (34) 1 (715) 
" ~2 miSu [2 4] k\-k, (7 - 15) 
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2. The (O g , 2+ , 3~ , 4 _ ) Configuration 



For this configuration, we simply quote the main results, detail, we now proceed by 
simply quoting the main results. The integrand evaluated on the solutions yields, 

X 6 =X 1 =X 2 = 0, X 5 =X 3 =X 4 = g3 4 (34)3 ( ^ ll|2] , (7.16) 

valid for any number of fermions and scalars. 

Because the integrands are at most linear in momenta the poles will be simple. It is then 
straightforward to evaluate the residues in each case by plugging in the leading momentum 
behavior at the pole. The only nonvanishing residues are 

Kg 4 - -Kg 2 - S 34 , K' A( ) 

or alternatively 

R = 4 (3 ^i) 32] (°» °> x > °> °> °> °) • ( 7 - 18 ) 

The integral coefficients are then given by 

r -P R 1 ,2 <34> 3 [32] 

and 

r - P R- V (34)3[32] 1 (720) 
° 2 ~ P * ' R ~ '2 S ^ (24) klh ■ (7 - 20) 

5. TTie (O g , 2" , 3+ , 4~ ) Configuration 

For this configuration, the results for the various solutions are, 

X 6 = Xx = X 2 = 0, X 3 =X 4 = [2|1^|3] [34] (42) 3 , 

r , lir 1; , 3 , ,(21^1 131 (31 £ 2 14) (42) (21 1^4 12) 
X 5 = [2)1413] [34] (42) 3 + (n f - 4) V 1 11 JV 1 21 M — ^ — 

•§34 

, 0(n n , 1 , (3|£ 2 |4) 2 (2|^ 1 |2)[2|1^ 2 |2) 

+2(n s -n/ + l) -2 (7.21) 

6 34 

where we have organized the X 5 expression such that first term is the Af = 4 sYM contribu- 
tion, the second term is the contribution from a chiral M = 1 multiplet, and the last term is 
the scalar contribution. Note that the fermions and scalars only propagate in the rightmost 
loop, so any Yukawa couplings are not relevant here. 
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All integrands but X5 are linear in momenta, giving simple poles. The X5 expression 
has up to fourth order poles. Somewhat surprisingly, upon careful analysis the residue of 
the simple poles receives no contribution from the fermions and scalars. All nonvanishing 
residues are again equal, 

Rg 4 = -Rg 2 = s u [23](A2) 3 [34\ , (7.22) 

or alternatively 

R = s u [2 3] (4 2) 3 [3 4] (0, -1, 0, 1, 0, 0, 0, 0) . (7.23) 
The integral coefficients are 

ci = Pi-i? = ^ 34 [2 3](4 2) 3 [34] , (7.24) 



and 



C2 = P 2 . R = _1 S34 [2 3] (4 2) 3 [3 4] ^L- . (7.25) 



4- The (O g ,2 ,3 ,4 + ) Configuration 

For this configuration the integrands are, 

X 6 = X 1 = X 2 = , X 5 = X 4 = [2|14|4] [43] (3 2) 3 , 
X 3 = [2|l£ 1 |4][43](3 2) 3 + (n / -4) 



3 , /A (2|£ 1 |4][4|£ 2 |3)(32)[2|1^ 2 |2) 



S34 

M9{ , 1 , [4|£ 2 |3) 2 (2|£ 2 £ 1 |2)[2|1£ 1 £ 2 |2) 

+2{n s — nf + 1) 2 • ('-^Dj 

s 34 

A helpful observation is that these integrands may be related directly to integrand results 
for the (O g , 2~, 3 + , 4~) configuration, by swapping k 3 A; 4 , X 3 -B- X 5 (and X x -H- X 6 ). 
This is a quite general feature as it follows from using the BCJ H —e nation 
(£2 - k 4 ) 2 A^(3, 4, -£ 2 , 4) = (£ 2 - fc 3 ) 2 A(°)(4, 3, -£ 2 , 4) for the product of trees that k 3 and 
k 4 attach to. 

As before the residues have no contribution from the fermions and scalars, and the non- 
vanishing residues are equal, 

p _ p _ ^ 34 (2 3) 3 [3 4] 

RQa - ~RQi - 77^\ > l'- 2 'J 

or alternatively, 

= ^2*||!M( 0i -1, 0, 1, 0, 0, 0, 0) . (7.28) 
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The integral coefficients are, 

r p R 1 323*34 (2 3) 3 [3 4] 

c 1 =P 1 -R=- ^4) ' (7 ' 29) 

and 

.i^ss^M i 

2 (24) fc>.jfe 4 K ' 

Finally, we note that the remaining nonvanishing helicity configurations (O g , 2 + , 3 + , 4 + ), 
{O g , 2~, 3 + , 4 + ), (C? ff , 2 + , 3~, 4 + ) and (0 9 , 2 + , 3 + , 4~) can be obtained by conjugating spinors 
in the results quoted above: that is using (ij) -b- [ji]. 



B. Scalar operator 

The scalar operator O s = ^Tr(0i 2 0i2) carries i?-charge. To obtain a non-zero result, 
the external states of the form factor must carry the opposite charge (in the all-outgoing 
convention). This means that the form factors involving only gluon states vanish. At three 
points only the form factor with two scalars, or one scalar and two fermions, are non- 
vanishing. Here we consider only the form factors for two external scalars and a gluon. The 
results are valid for any massless theory with rif > 1 fermions and n s > 1 scalars that has 
a Yukawa coupling gy u kawa'4' 1 i j2( Pi2, and (/Yukawa = 9ym- Then all the cuts become identical 
to those of M = 4 SYM. 

1. The {O s , 2+ ,3+, 4+) Configuration 

For the configuration (O s , 2+, 3+, 4 + ) the cut integrands are 

X 6 = X l = X 2 = , J 5 = J 3 = J 4 = (3|(£x + fc 4 )|2] (2 3) [3 4] 2 (7.31) 

All the poles will be at most simple, and it is straightforward to obtain the residues in each 
case. The only nonvanishing ones are 

7? - 7? -s s [M] (23) (7V) 

Kg 4 - -Kg 2 - S23S34 , { I -oZ) 



or, using vector notation, 



R = ^23^34 ^^ (0, -1, 0, 1, 0, 0, 0, 0) . (7.33) 
31 



The integral coefficients are then given by 



and 



c 1 = P 1 -fl=i» a »«M|5>, (7.34) 



2. The (O s ,2+ ,3+, 4+) Configuration 

The case (C s , 2+, 3 + , 4+) is related to above configuration via the swap fc 3 A; 4 , X 3 -h- X 5 , 
Xi -H- X 6 . The integrand is then 

X 6 = Xi = X 2 = , X 5 = X 3 = X 4 = (4| (£x + fc 3 ) |2] (2 4) [4 3] 2 (7.36) 

The nonvanishing residues are again equal 

d _ d _ „ „ [34] (2 4) 
itft, - --Kg 2 - -S23S34 — — , 



alternatively we have 



R = -^ J 3 ^ 4 ^ , -1, 0, 1, 0, 0, 0, 0) . (7.38) 



The integral coefficients are then given by 



and 



c 1 = P 1 -R= -I a2sas4 2||^ , (7.39) 



5. The (O s ,2 + ,3f,4f) Configuration 



Finally, for the configuration (O s , 2 + , 3+, 4+), we have 



2^=^ = 0, i 5 .x 3 ^x^Wl+|M£li. (7 . 41) 



And the nonvanishing residues are 



, [2 3] (3 4) 



Rg 4 = -Rg 2 = s 34 — jp-^j — , (7.42) 
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alternatively 



, [2 3] <3 4>_ 



R = ^4 ' (0, -1, 0, 1, 0, 0, 0, 0) . (7.43) 



and 



The integral coefficients are then given by 

^fl.^-^MM 1 . (7.45) 
2 34 (24) j^.fc 4 v ' 

The remaining nonvanishing gluon-scalar configurations (O s , 2+, 3+, 4~), (C s , 2+, 3~, 4+) 
and (O s , 2~, 3+ , 4+) can be obtained from spinor conjugations of the above ones. (The fact 
that the scalars are complex does not interfere with this relation, as has been explicitly 
checked.) 



4- Comparison to Previous Results for (O s , 2f, 3f, 4 + ) 

Here we compare to the results of Brandhuber, Travaglini and Yang [79|. They obtain 
the following answer for the numerator of the double box integral: 

The integrand of the maximal cut is simply this numerator Zj = N^-q for all kinematical 
solutions. Clearly this integrand is different from the ones we computed above. However, 
the difference should be due to terms that integrate to zero. 

Computing the residues we recognize the same expression that appeared in eq. (17.32p . 
But this time there are more nonzero residues, 

R - R -R - R - J < , t 34 J < 23 > (7A7) 
Kg 3 - —tig 1 — Hg i — -Kg 2 - — S 23 S34 — J^j — • {' A ') 

Two of the residues vanish due to an identity, 

Rg 5 = Rg 6 = ^ 3 (2 4) 3) (2s24fc i ' h ~ 2$2 ^ ' h) = ' (7 ' 48) 
Writing this as a vector gives 

R = ^s 23 s J 3A ^ 3) (-l, -1, 1, 1, 0, 0, 0, 0) . (7.49) 
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The integral coefficients are then given by 

Cl = p i .R = 

and 



1 [34] (2 3) 

-S23>S34- 



C 2 



P2 ■ R — — -S23S34 



(2 4) 
[3 4] (2 3) 



(7.50) 



(7.51) 



(24) k\-k 4 ' 

These results agree with the integral coefficients that we computed in section I VII B 1[ so 
we find complete agreement with the results of ref. 79]. 



VIII. CONCLUSIONS 



In this paper, we have continued the approach of ref. [61] in developing the maximal 
generalized unitarity method for two-loop amplitudes. Cutting propagators can be viewed 
as replacing the original real loop-momentum contours of integration by contours encircling 
the global poles of the integrand. Double-box integrals, independent of the configuration of 
external masses, have eight independent global poles. The possible ways of cutting propa- 
gators, or equivalently the allowed linear combinations of contours surrounding global poles, 
are constrained by the requirement that the evaluation along any contour respect the van- 
ishing of certain insertions of Levi-Civita symbols, as well as of total derivatives arising from 
integration-by-parts identities. We derived the corresponding constraint equations for dou- 
ble boxes with up to three external masses. The number of master integrals depends on the 
configuration of the external momenta, but in all cases the constraint equation yield unique 
and simple formulae for the coefficients of the master integrals. The master formula (16. ip . 
along with these projectors, given in eqs. ( 16 . 3H and (I6.5p . are our principal results. 
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